Self-accelerating Massive Gravity: 
Exact solutions for any isotropic matter distribution 
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We present an exact solution to the equations of massive gravity that displays cosmological 
constant-like behavior for any spherically symmetric distribution of matter, including arbitrary 
time dependence. On this solution, the new degrees of freedom from the massive graviton generate 
a cosmological constant-like contribution to stress-energy that does not interact directly with other 
matter sources. When the effective cosmological constant contribution dominates over other sources 
of stress energy the cosmological expansion self-accelerates, even when no other dark-energy-like 
ingredients are present. The new degrees of freedom introduced by giving the graviton the mass do 
not respond to arbitrarily large radial or homogeneous perturbations from other matter fields on 
this solution. We comment on possible implications of this result. 



I. INTRODUCTION 

More than seventy years have elapsed since Pauli and 
Fierz made the first attempt at writing a theory of grav- 
ity with a massive graviton In the intervening years, 
daunting challenges to realizing such a theory have been 
found, including the scylla of incompatibility with Solar 
System tests 0, [31 and the charybdis of ghost- like de- 
grees of freedom Q. Recently, de Rham, Gabadadzc, 
and Tolley have constructed a theory of massive grav- 
ity [EHl] that evades these dangers [9, [HI ■ This theory 
also contains a vacuum solution that recovers exactly a 
Schwarzschild-de Sitter solution [ll|, E3- Moreover, in 
the flat matter dominated limit, the theory has a solution 
that responds to the presence of matter by producing an 
effective cosmological constant contribution to the stress 
tensor at the cost of introducing inhomogeneous solutions 
for the Stiickelberg fields that describe the new degrees of 
freedom that come from massive gravity ff3T | . For an open 
universe, a related solution has been explicitly shown to 
evolve into self-acceleration fl4| . 

In this paper, we generalize considerations in fl3l Ibij 
to an arbitrary spatially isotropic metric. We find cos- 
mological constant type solutions in the presence of any 
isotropic distribution of matter. Such solutions connect 
the fiat matter dominated solution [l3| to the de Sitter 
solution [Tl|, [l2| allowing a cosmological expansion his- 
tory identical to the ACDM model even in the presence 
of spherically symmetric matter perturbations. 



II. MASSIVE GRAVITY 

The covariant Lagrangian density for a theory of mas- 
sive gravity will have, in addition to the usual Einstein- 



Hilbert term, a mass term with a potential U, 



in' , 



(1) 



M p 2 , = 1 



K^v is a ten- 



ttG and H = c = 1 throughout, 
sor that characterizes metric fluctuations away from a 
fiducial (flat) space time. At the linearized level, the po- 
tential must take on the Fierz-Pauli structure to be ghost 
free; but any purely linear theory exhibits the vDVZ dis- 
continuity [2j, |3| , where an extra helicity mode couples to 
matter even in the m — > limit. Nonlinear extensions 
to the Fierz-Pauli potential can evade this problem via 
a strong coupling phenomenon known as the Vainshtein 
mechanism [151 ] . where the extra coupling is suppressed 
near matter sources. However, these extensions typically 
contain an unhealthy ghost-like degree of freedom Q. 

For a theory of massive gravity to be free from this 
ghost, the potential term must take a special form built 
out of expressions that have the form of total derivatives 
in absence of dynamics [1] . These can be written as con- 
tractions of the tensor 
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The matrix VS is understood to denote 

The potential-generating matrix is defined as 

E"„ = g^d a rd^ b Vab = g^ av , (3) 

where <f> a are the 4 Stiickelberg fields introduced to re- 
store diffeomorphism invariance. The 4> a fields transform 
as scalars, while S, and K. transform as tensors under 
general coordinate transforms. 

In matrix notation, the potential can be written 0, Q 

-U/4 = [K] 2 - [K 2 ] + a 3 ([/C] 3 - i[K][K 2 ] + 2[7C 3 ]) 
+ a 4 ([K} 4 - 6[7C] 2 [7C 2 ] + 8[7C][7C 3 ] + 3[7C 2 ] 2 
-6[7C 4 ]), (4) 

where brackets denote traces, [A] = A^^, and a^, a 4 
are free parameters. Using Eq. @, we can reexpress the 
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potential in terms of traces of products of v E 
j = -12 + 6[V£] + [£]-[\/E] 2 

+ a 3 (-24 + 18[V£] - 6[\/S] 2 + [V£] 3 

-3[E]([7S]-2) + 2[S 3 / 2 ]) 
+ a 4 (-24 + 24[VS] - 12[V£] 2 - 12[\/£][£] 
+ 6[%/S] 2 [£] + 4[/S] 3 + 12[£] - 3[E] 2 
-8[E 3 / 2 ]([V£]-l) + 6[£ 2 ]-[\/£] 4 Y (5) 



Variation of the action with respect to the metric yields 
the modified Einstein equations 



1 



G — m 2 ^^ A 7^™) 



(6) 
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where G^ is the usual Einstein tensor and T^' is the 
matter stress energy tensor. Here 

— U 



T (K.) = _i 5 _, 

= - \{^9,u - 2£ M „ - 2(3 - 
+ a 3 \ - 3 (6 - 4[VE] + [%/S] 2 - [S] 



(7) 



6 f [VS] - 2) - 6E 3 / 2 



a 4 [-24 (s 2 „ - ([VS] - 1)E 3 {, 2 ) 

- 12 (2 - 2[VE] - [53] + [\/E] 2 ) £ M „ 

- (24 - 24[VE] + 12[\/E] 2 - 4[v / E] 3 

-12[53] + 12[53][\/£] - 8[E 3 / 2 ]) Ve„„] } 

is the dimensionless effective stress energy tensor pro- 
vided by the mass term. Note that this effective stress 
energy depends explicitly on the metric itself. To solve 
the modified Einstein equation, we first parameterize the 
metric and then solve for the joint effect of the matter 
and mass term. 

III. EXACT SOLUTION 

Generalizing [H,[ll[, we consider an arbitrary spatially 
isotropic metric, 



ds 2 



b 2 (r, t)dt 2 + a 2 (r, t)(dr 2 + r 2 dVt 2 ). (8) 



We correspondingly take a spherically symmetric ansatz 
for the Stiickelberg fields: 



<t>° = f(t,r), 



4?=g{t,r)-, 
r 



(9) 



and look for solutions to the functions g(t, r) and f(t, r). 
The potential matrix ([3]) then takes the form 



£ = 



//' - 99' 



b 2 b 2 

gg'-ff -f 2 + g' 2 
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(10) 



where primes denote derivatives with respect to r and 
overdots with respect to t. 

The resulting, rather involved, calculation is made eas- 
ier by isolating the upper-left-hand 2x2 submatrix of 
E and using the Caylcy-Hamilton theorem, which states 
that a matrix solves its own characteristic polynomial. 
For a 2 x 2 matrix A, this means 

[A] A = A 2 + (dct A) I 2 , 

where I2 is the 2x2 identity matrix. We can then use 
that detA™ = (detA)™ to find the square root of £2, 
the upper- left- hand 2x2 submatrix of E: 
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where 



X 
W 



'\ 2 
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(11) 
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(12) 



and n = sgn(/.g' - gf). 

With Eq. JlT]), traces of E" become 



E] = 
[E] = X - 2W 



2g 

ar ' 



(13) 



2g 2 



[E 3/2 ] =X 3/2 -3WVX 



2g 3 



[53 2 ] = X 2 - 2W(2X - W) 
and the potential is given by 

^ = p (JL) + Vx Pl (^ 

4 \ar/ \ar/ 

where the P n polynomials are 



WP 2 (±) , (14) 



P (x) = -12 - 2x(x - 6) - 12(x - l){x - 2)a 3 

- 24(x - l) 2 a 4 , 
Pi(x) = 2(3 - 2x) + 6(x - l)(x - 3)a 3 + 24(a; - l) 2 a 4 , 
P 2 {x) = -2 + 12(2! - l)a 3 - 2A{x - l) 2 a 4 . (15) 
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Varying the action with respect to / and g yields the 
Stiickelberg field equations 



a 3 r 2 



I 

a 2 br 2 



P 1 + p,a 2 r 2 g'P 2 



(16) 



9 . ?\v 



and 



a 3 r 2 
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X V6 
a 2 br 2 



fj,a 2 r 2 gP 2 



2 2//n 



^ + -)Pi+^a 2 r 2 fP 2 
b a / 



z 2 6r 



(17) 



where Pn(x) = dP n /dx = ardP/dg. By inspection, we 
find that a solution to the / equation of motion, Eq. (|16|) , 
is given by Pi(xq) = 0, or 



separately tracking the equal 9 and <j) contributions to 
2(g/ar) n terms in the traces of Eq. (TIB")) . Their separate 
variations can then be reduced with Eq. (TH))) . 

Hence, the effective energy density and pressure are 



(m 2 Ml x )T^ K) = 



-Pic 
vk 
p K 
PK 



where 



Pic = ~PK = 2 m2A// pi- p o(^o)- 



(22) 



(23) 



This shows that a cosmological constant type solution 
exists for general isotropic metrics. Conversely, the mod- 
ified Einstein equation for arbitrary spherically symmet- 
ric distributions of matter becomes the ordinary Einstein 
equation plus a cosmological constant on this solution. 

For example, the spatially flat FRW space-time is a 
subset where a(r, t) = a{t) is the scale factor, b{r, t) = 1 
and the modified Einstein equation (J6j) just becomes the 
usual Friedmann equation 



1 + 6a 3 + 12a 4 ± yj\ + 3a 3 + 9a 2 - 12a 4 

3(«3 + 4d!4) 



(18) 



and hence g = xoar. Note that if 03 = ay = 0, -Pi (2;) 
becomes linear and g = Zar/2 is the solution. 

The equation of motion for g evaluated on the solution 
provides a constraint on / 



XP[ 



m 

x 



p^w-p^, 



(19) 



where the P n functions are evaluated at xq and we have 
used the fact that 



W 



(f+-rf--rf')x . 
V a a / 



(20) 



An explicit solution for / is not required for the com- 
putation of the stress energy tensor. That is, the physical 
background solution does not depend on the choice of so- 
lution for / and in particular is independent of the spatial 
and temporal integration constants that are introduced 
in solving for /. After straightforward but tedious alge- 
bra, we find that its nonzero components are: 

T& = lPo(x )b 2 , 
T% = -^P (x )a 2 , 



rp A 

J- a, 



2^2 



-P (x )a 2 r 



(21) 



The Tq5 and T^ r pieces can be easily checked from 
Eq. p^|) by direct variation with respect to g tt and g rr , 
noting that the polynomial pieces come from the angular 
metric. The angular pieces can be similarly analyzed by 
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Mil, 



{PK + Prn)- 



(24) 



FRW space-times with spatial curvature K ^ are also 
included with 



a(r,t) -> 



i(t) 



(25) 



in isotropic coordinates. Note that for the FRW metric, 
this solution applies for radiation and matter domination 
as well as for a self-accelerated epoch where the mas- 
sive graviton itself provides the cosmological constant- 
like dark energy. It also allows for arbitrary isotropic 
perturbations around the FRW metric with 



a 2 {r,t) =a 2 (t)[l + 2$(r,t)], 
6 2 M)=[l + 2vI/(r,t)]. 



(26) 



Thus the solution remains of the cosmological constant 
type for arbitrary spherically symmetric matter distribu- 
tions. Furthermore, the matter only sees the effects of 
the mass term as a cosmological constant with no direct 
coupling to the Stiickelberg fields on the exact solution. 

It is straightforward to verify that our class of solu- 
tions subsumes several particular solutions that have pre- 
viously appeared in the literature: in vacuum (#=0), 
it recovers exactly the static Schwarzschild-de Sitter so- 
lution from [HI, [IH an d fl6l : it also reproduces the decou- 
pling limit solution in Il7l . ll8| [ as well as the open universe 
solution reported in |l4j . Reduction to these particular 
solutions is effectively made through a choice of f(t,0), 
thus demonstrating that those solutions are not unique. 
Our solution is also similar to another Schwarzschild-de 
Sitter solution [l9| . as well as the vierbein formulation 
solution of 1201. 
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IV. DISCUSSION 

The solution we have found is a perfect analog for a 
cosmological constant. Because the solution exists for 
any isotropic distribution of matter, it recovers static so- 
lutions like Schwarzschild-de Sitter in vacuum and gener- 
alizes them to dynamical cases such as the FRW cosmol- 
ogy. In each of these cases, the presence of other isotropic 
sources of stress-energy does not alter the cosmologi- 
cal constant-like behavior of massive gravity. Hence, 
we can have a truly self-accelerating gravitational back- 
ground that coexists peacefully with a standard cosmo- 
logical history; the self acceleration begins in precisely 
the same manner as cosmological-constant-driven accel- 
eration would begin, only here the size of the apparent 
cosmological constant is set by the graviton mass and the 
other free parameters of the theory (0:3 and a.±). More- 
over, the extra gravitational degrees of freedom present 
in this theory do not appear to couple to radial matter 
perturbations. In practice, this means that radial matter 
perturbations will feel only the ordinary gravitational at- 
traction as in general relativity. This is in contrast with 
the enhanced gravitational force felt by matter perturba- 
tions around other solutions of this theory. Likewise the 
Stiickelberg-driven self-accelerating background physics 
will not respond to spherically symmetric matter pertur- 
bations either. 

In this paper, we have restricted ourselves to isotropic 
situations. Note that although we have assumed our 
Stiickelberg fields to be in a radially symmetric configura- 
tion, their effective center in space disappears in their ef- 
fective stress energy, which is homogeneous and isotropic. 
Indeed, we can recover fully homogeneous background 
solutions supported by the Stiickelberg fields. This sug- 



gests that perhaps even more general inhomogeneity in 
the matter fields may not drive inhomogeneity in the ob- 
servable effective stress energy of the Stiickelberg fields. 

The obvious next step in assessing the solution we have 
found is to attempt to study perturbations around it. 
Naively speaking, the solution we have found gives us 
reason for concern, because the part of the action that we 
might expect to generate kinetic terms for these fluctua- 
tions appears to vanish on our solution, since Pi(a;o) = 0. 
At first, this appears to confirm the findings of plj . 
who find vanishing kinetic terms for perturbations of the 
Stiickelberg fields around their open Universe solution 
14|, which is a member of our class of solutions. How- 
ever, our class also includes the solutions found in fl8l ]. 
who studied vector and scalar perturbations in the de- 
coupling limit of their solution and found non-vanishing 
kinetic terms for the scalar perturbations as well as ei- 
ther ghost-like or vanishing kinetic terms for the vector 
perturbations. Similarly, the closely related decoupling 
limit solution of [I?) also has non- vanishing kinetic terms 
for scalar perturbations. In light of these considerations, 
a careful study of general perturbations to our class of 
solutions will be an important area for future work. 
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